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Abstract 

The off-diagonal heat-kernel expansion of a Laplace operator including a general gauge-connection 
! is computed on a compact manifold without boundary up to third order in the curvatures. These 

lO . results are used to study the early-time expansion of the traced heat-kernel on the space of 



OO 



transverse vector fields satisfying the differential constraint D^v^ = 0. It is shown that the 
. resulting Seeley-deWitt coefficients generically develop singularities, which vanish if the metric is 

flat or satisfies the Einstein condition. The implications of our findings for the evaluation of the 
gravitational functional renormalization group equation are briefiy discussed. 
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1 Introduction 



The heat-kernel technique is a mathematical tool with a wide range of applicability both in 
mathematics and theoretical physics [1, 2, 3]. This is owed to the fact that it provides a formal 
way to treat functional traces and determinants of local differential operators, including any 
Laplace-type (second order differential) operator. Irrespective of the particular context they 
appear in, the arguments of such traces can often be written in an exponentiated form, that 
allows their systematic treatment in terms of the heat-kernel 

H{s) = e-^^o . (1.1) 

Here Aq denotes a generalized Laplace operator on a closed and torsionless Riemannian manifold 
that we will refer to as "spacetime". In the context of theoretical physics, it allows to give, for 
example, a definition of the propagator^ and 1-loop effective actions.^ Moreover, it allows to 
compute counterterms and anomalies in an elegant way that can be extended naturally to field 
theories on arbitrary curved space. Furthermore, it constitutes an essential ingredient in solving 
the gravitational functional renormalization group equation [4], see [5, 6] for reviews. 

The trace of the heat-kernel has two main expansion schemes. One is the local or early-time 
expansion in terms of powers of s 

1^^(s) = — / d^x^s^tiA:., (1.2) 



n>0 



generally referred to as Seeley-deWitt expansion. The scalar quantities An{x) are called heat- 
kernel coefficients and are local functions of the curvature invariants and their covariant deriva- 
tives. The other scheme is based on a non-local expansion in terms of curvature tensors and 
schematically represented by [7] 

1 r " 

Tr H{s) = —-^ E ^ / n (^'^.V^) F (sA,„ . . . , 7^,, . . . 7^,„ . (1.3) 

This expansion involves arbitrary powers of derivatives at every order in the curvature tensors 
TZi = TZi{xi). In (1.3) each operator Aj is acting only on the corresponding invariant TZi and the 
tensor structure of the invariants has been suppressed for brevity. Throughout this paper, we 
will work with the early-time expansion (1.2), since its coefficients and their derivatives can be 
computed recursively. 

In quantum field theory (QFT), the computation of (1.1) plays a central role, since it allows 
explicitly covariant computations of which the aforementioned propagator and 1-loop effective 
action are only special examples [8]. In fact, a general trace involving a function of Aq can 
formally be related to the traced heat-kernel by a Laplace transform 

/•oo 

Tr/(Ao)= / dsfis) TrHis). (1.4) 
Jo 



^ Using the well known Schwinger-trick — ds e 
^ri.ioop = iTrlog A = i /J- i^Trfe-'^l. 
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Thus, a wide class of covariant computations in QFT, involving different functions /(x) of the 
same differential operator Aq, are reduced to the calculation of the single object Tr H{s). 

The caveat of this approach is that the typical operators Aq of physical interest are often 
very complicated and, in particular, not of Laplace-type. A prototypical example is Yang-Mills 
theory where the gauge fixed inverse propagator of a gauge field A'^^ is given by a differential 
operator that has indices on the internal gauge group as well as spacetime, and, for a general 
gauge-parameter a, is of non-minimal form in the sense that it includes uncontracted derivatives 

A^^. = -g.^D^'D'^"' + (1 - a-') D^dI' , with Df = 8^6'^'' + r^'^A'^ . (1.5) 

One possibility to simplify such types of operators is the use of covariant projectors [9, 10, 11] 
which leads to Laplace-type differential operators acting on subspaces of the original field space. 
A typical example for such projections are the transverse decomposition of a vector field, or 
the transverse traceless decomposition of a symmetric tensor field [12]. Alternatively, one may 
carry out a resummation of the non-minimal derivative terms in the operator [2, 13]. These 
two methods share a common building-block, the generalized heat-kernel traces which contain 
insertions of covariant derivatives 

(1.6) 

Traces of this type have been studied in [11, 14, 15, 16] and their Seeley-deWitt coefficients can 
be obtained by a generalization of the early-time expansion (1.2). In particular, they prove very 
valuable when evaluating traces of operators where the heat-kernel expansion (1.1) is not known. 

The rest of the paper is structured as follows. In section 2 we study traces of the form (1.6), 
using the deWitt-algorithm for determining their curvature expansion recursively. Implementing 
this algorithm in a computer algebra system [17, 18], we generalize the results of [11, 14, 15, 16] 
to differential operators on a general gauge bundle including an arbitrary endomorphism. As an 
application of these results, we study the early-time expansion of the traced heat-kernel of the 
Laplace operator on the space of transversal vector fields in section 3, giving the Seeley-deWitt 
coefficients up to order 7i? in table 2. We close by briefiy commenting on possible applications 
of our findings in the context of QFT in section 4. Some technical details have been relegated to 
four appendices. 



Tr 



H{s) 



2 The ofF-diagonal heat-kernel 

In this section, we outline how the derivatives of the heat-kernel coefficients (1.2) can be computed 
recursively. The main virtue of the off-diagonal method used here is that it allows to compute 
operator traces of the general non-minimal form (1.6), which have a wide range of applications 
in QFT. 



3 



2.1 Recurrence relation for off-diagonal heat-kernel coefficients 

We assume that our spacetime is a closed Riemannian manifold without boundary and of arbitrary 
dimension d. The Laplace operator Aq in the heat-kernel (1.1) is taken to be of general second 
order form 



^0 = - g^'d^d, + a^d^ + b 



(2.1) 



where in the second line it is cast into standard notation [1], involving a covariant derivative 
operator = + and an endomorphism E.^ The symbol A is reserved for the Laplacian 
built from the covariant derivative only A = —D^^D^ (without an endomorphism) and is 
the covariant torsionless spacetime derivative compatible with the metric g^^- We define to 
be an unspecified connection on an internal bundle over the spacetime manifold and we assume, 
without loss of generality, that E is an endomorphism on the same bundle of A^] whenever this 
is not the case it is sufficient to decompose their bundle into the direct sum of their respective 
bundles. Further we define the sum of the curvatures of the connection A^ and the Levi-Civita 
connection 

F^,u = d^A^ - d^Af, + [A^, A^] + [V^, V^] . (2.2) 

For notational simplicity, all internal indices are suppressed, so that the quantities A^, F^^j and 
E are understood as matrices on the internal space. 

The heat-kernel H{s) owes its name to the fact that it fulfills a generalized heat-equation 
with boundary condition i^(0) = 1, where 1 denotes the identity on the internal space. To arrive 
at the heat-equation in the form of an ordinary partial differential equation, it is convenient to 
express H{s) in terms of its matrix elements, called the off-diagonal heat-kernel 

H{x,y;s) = {y\H{s)\x) = (?/|e-^^°|x) , (2.3) 

since the basis elements are at different points x and y of the manifold. The definition (2.3) is 
equivalently given as an initial value problem of the heat-equation 

(9, + Ao,..)^^(x,y;s) =0, H{x,y-Q)=5^^yl. (2.4) 

The solution of this differential equation has the interpretation of heat propagating on the man- 
ifold, according to the operator Ao,xi from a source located in y. Here s plays the role of a 
diffusion-time for the process. 

The initial value problem (2.4) can be solved explicitely in the simple case of a flat manifold, 
where both the connection A^ and the endomorphism E vanish. In this case the solution of (2.4) 
is given by 



1 _ (x-yr 
(47rs)' 



H{x,y;s)\^,, = ^^—^e 4s . (2.5) 



''This can always be done as long as the manifold is torsionless [19]. 



4 



As an ansatz for the general solution, this expression is modified by introducing a general function 
n{x,y;s) 

H{x,y;s) = -^— -e-^ nix, y;s). (2.6) 
(Atts) ' 

Here a{x,y) is half of the squared geodesic distance, called the "world function" [8]. It generalizes 
the distance measure in (2.5) in a covariant way and satisfies 

^ V''^ = ^ (2-7) 

for arbitrary spacetime points x and y.^ In order to find the partial differential equation satisfied 
by n{x,y;s), one substitutes the ansatz (2.6) into the heat-equation (2.4) to obtain 

ids + Ao,.) H = -J—e^i (--^n + + -(T;^f)/' + O^n - J];/ + En) . (2.8) 

[Airs) ' \ zs s J 

The heat-equation is solved if the bracket on the right hand side vanishes. Inserting the early-time 
expansion of the heat-kernel 

n{x,y;s) = Y,s^An{x,y), (2.9) 

n>0 

and requiring that the bracket in (2.8) vanishes at all orders of s independently, this equation 
yields the master-equation 

(n - I + icT;/) A„ + a-'^A^;^ - + = , n>0, (2.10) 

subject to the initial conditions = and = 1- 

Equation (2.10) still constitutes a complicate partial differential equation for the coefficients 
An{x,y). In order to solve it, we exploit that the full off-diagonal heat-kernel coefficients for 
non-coinciding points can be expressed as the geodesic expansion [14] 

An{x, y) = Yl M- ^/^i • • • ^M-^" • • • ' (2.11) 
z — / m\ 

where the overline denotes the coincidence limit of any bi-tensor C{x,y) 



C{x,y) = C{x,x). (2.12) 



Notice that here, it is sufficient to know the quantities D(^^^ . . . D^^^An symmetrized in their 
indices.^ Substituting (2.11) into (2.10) then allows to recursively determine the expansion- 
coefficients (2.11) at coincident points y — t- x. Since only covariant expressions are used in 
its derivation, the An{x,y) are given in terms of an expansion in curvature monomials. To 
solve (2.10) for any An, the coincidence limits of the derivatives of and a are required. 

Comparing powers of the curvatures occurring in these objects, the systematics is easily found 
and summarized in table 1. For example, to compute all coefficients up to the 6-derivative order 
(7^^), one needs the coincidence limit of 8 derivatives acting on a, 6 derivatives acting on Aq and 
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Table 1: Analysis of the terms entering the recursion relation (2.10) at coincidence limit (we omit 
the overline for brevity). In order to compute an entry, one needs to compute every object that 
is above and on its left. Here TZ counts the number of curvature tensors and of square covariant 
derivatives. (For example D^R^ counts 7^"+*"/^.) 

so on. The (t-^^,,,^^ are found with the help of the defining equation (2.7) [20]. With the initial 
condition a = 0, it is straightforward to find the first few expressions, by successively applying 
derivatives to the equation. The coincidence limits up to fifth order in the derivatives are 

and in addition we give 

-= l^I? 4p 7e> 4c> E> 7_| 8_r) jy 7 

^■,aP{pvpa) ■^{p\a\u\l3\;pcr) ^^{^u\a\ ^p|/3|o-)7 ^^'){^u\a\-^pcr)l3 ~r i'!^-^^{pu\a\^p\l3\u) 

_l_16r) p 7 8.f?, 7_i_4p f?,7 

' ii^-'^^{pup^a)aj3 \^-'^'y[pup^a)f}a ' i^^"j{p,vp^a) 



(2.14) 



ah 



These objects can be found as follows. After taking n derivatives of (2.7), terms appearing with 
n + 1 derivatives vanish in the coincidence limit because does. Inserting all lower order results 
and commuting the n indices in the remaining terms into a unique order immediately reveals 
the result. The method is conceptually straightforward and requires only simple algebra, yet the 
higher order terms increase in size quickly and render a by-hand-computation unfeasible. 



Equation (2.10) is solved for the quantities D^^ . . . D^^An for any n, m > in the same way. 
One can obtain an infinite set of algebraic equations for these by applying derivatives and taking 
the coincidence limit of (2.10). These can be solved for recursively, substituting all lower order 
objects to find the next order. Once all required derivatives are found for some n, one can proceed 
to n + 1 until all ingredients to the heat-kernel expansion to a desired order are found. The results 
are listed in the next subsection. 

Notably, the recursive relation for the heat-kernel coefficients (2.10) becomes independent of 



the spacetime dimension once the coincidence limit is taken, since cJ;^'^ = d cancels the multi- 
plicative d in (2.10). This is an important observation because it implies that the heat-kernel 



We use a semicolon to abbreviate any covariant derivative, D^^a = a.^. We use the convention that all 
derivatives are with respect to the coordinate x. 

^It is important to note that in general -D(pi . . . -D^^^jAn ^ • • ■ ^Min)^"- 
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coefficients of Laplace-type operators cannot depend on the dimension explicitly. We stress that 
this property does not hold for more general operators like (1.5). 



2.2 Heat-kernel coefficients on a general vector bundle 

Starting from the discussion of the previous subsection, it is straightforward to implement the 
recursive equations (2.7) and (2.10) in a computer algebra software [17, 18] to find the ^n;^i..^tm 
explicitly. Up to second order in the curvatures, the coincidence limit of the coefficients and their 
derivatives are found as 

= 1, 



^M^o = 0, 
1 



3 1 1 

A = -e + \r, 

'd;^! = -\E,^.-\F,^.■:' + ^R;^., (2.15) 

L>(^Z)^)Ai = --E.[^^) - -R^uE - -F^{^,"u) + ■7Ea{uE°'f,) 



1 1 1 

20 '^^'^> 60 ^ 36 



30 72 180 ^ 180 

We also computed all coefficients up to third order in the curvatures (7^"^). Here we state only 
the explicit result 

= -1^3 _ ]_f^E)E + —E.,,E.'' - —E(AE) - —(AAE) 
d 6 12^ ^ 12 ' 12 ^ ^ 60^ ^ 

1 111 

J 7?M ppi^ I ppf p P J ppp p 1^ p p>^P pp 

180 "'^ ''' 45 ^ '^'^ 90 '^^ 90 ^ 

+—E^R + —(AE)R - —E.^R.^ + —E(AR) - —E-.^R'"' 
12 36^ ^ 30 • 30 ^ ^ 90 

-—ER^ + —ER.^Rf"' - —ERaupaR^"'"' (2.16) 
72 180 ^ 180 ^ ' 

_|_J_77 pP^r}J^}_p pPP-Q^ p T)P^P<^ 

72 30 180 
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+—(AAR) - —R(AR) + -^R^Rt" + —R.^^R'"' + —R^JAR^'n 
280^ ^ 280 ^ ' 5040 ' 420 630 ' 

Rpv;pR^ — -ir,nn.Rpy;pR '''^ ~ 77^ -^/i'^pcr ( Ai?^ ^ ) + -r^R,_iupa;\R^ ^ 



2520 ^"^"^ 1260 '^'^"^ 420 '^'^'^^ ^ ' 560 

1 3 1 1 

— TT^RRpuR^'^ + ccvn RpuR^^R'^ 



1296 1080 5670 

_ ^ p p Tjppucr _i_ ^ p p Tjpvpcr _ ^ p p pv 

1890 ^'^ '^ 1080 '^'"^ g^g-"'Mi'-"'WO"^-"' PO"^ 

J 3_ m P a R^"'^'^ -I RM*^ /? a RP'^'^I^ 

567 ^ o--n,^op/3-n, 11340 P(^^^^vap^ 

Some additional coefficients at third order in the curvatures are given in appendix C. The ones 
that are not given are too large to be conveniently written. This formula generalizes previous 
computations of A3 with traced internal space indices [1], and in the limit of trivial connection 
{F^y = 0) [14]. Note that in the above expressions, a full covariant derivative D is identical to 
the covariant spacetime derivative V when it acts on the spacetime curvature tensors R. 

2.3 Evaluating traces with non-minimal differential operator insertions 

The derivatives of the heat-kernel coefficients can be used to evaluate operator traces with non- 
minimal derivative insertions. The trace of the heat-kernel defined in terms of its matrix elements 
(2.3) is 



Tr[e-"^°] =try" d'^x^ {. 



x\e-'^° \x) 



tr J d'^x^ H{x,x;s) , 



(2.17) 



where "tr" is the trace in the internal space. This expression generalizes to the case of arbitrary 
insertions of uncontracted covariant derivatives like 



Tr[D^, . . . Z)^„e-^^o] =tr ^ d'^x^ {x\ D^, . . . D^^e 
= ii j d'^x^ H^^,„f,„{x,s) , 



(2.18) 



where we introduced the abbreviations for the coincidence limits 

Hp,...pAx, s) = D^,... D^^H{x, y- s) = {x\ D^,... I^^^e"^^" |x) . (2.19) 

The derivatives acting on H{x,y;s) can be applied explicitly to the definition (2.6), in order to 
express i?^^...^^ (x, s) in terms of the heat-kernel coefficients (2.15). Hereby the derivatives will 
act on both cr{x,y) and Q{x,y;s), where the latter is given by the heat-kernel coefficients via 
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(2.9). We give the result for the first six of these matrix elements with symmetrized derivatives: 

n>0 

H^{x, s) = (47rs)-'^/2 s^^;^^ 



n.>0 



n>0 



^(M^p)(a:,s) =(47rs) '^/^J2^''~H-^9{puD^)An + D^pD^D^)An-i 



n>0 



ri>0 



+ D(^xDpD^Dp^An-2 
15 



(2.20) 



n>0 



^9{axDpD^D^)An-i + D^^DxDf,D^,D^)An-2 
15 



H[^ypXaP) {X, S) = (47rs) '^/^ ^ s" ^ ^9{fia9Xp9i^p)An 



n>0 

45 

+ -r9{l3a9XpDuD^)An-l - —gi^i3aD\DpDj,D^-^An-2 



+ D(^pDaDxDpD^D^)An^3^ 

In these expressions, the coincidence limits of derivatives of a (2.13) have been inserted, and 
the sums are understood with the boundary conditions A_i = A^2 = • • • = 0. The general 
unsymmetrized formulas can always be recovered from the symmetrized ones by commuting the 
derivatives. For example we have 

Hp,{x,s) =H(^^,){x,s) + (47rs)-'^/2^s'^-iZ)[^Z),]A„„i 

, (2.21) 

=i/(^,) {x, s) + (47rs)-''/2 Yl s'^-'Fp.An-i , 

n>0 

where the commutator in the second term becomes the curvature tensor. With the explicit 
knowledge of the expansion coefficients (2.15), all traces of the form (2.18) are given. 

The relevance of these traces lies in the fact that quite general expressions involving the 
Laplacian operator and uncontracted derivatives can be computed by use of (2.20). For example 

Tr[L>^, . . . Dp J (Ao)] = tr y d'^x^ J ds f (s) D^, . . . Dp^H{x,y; s) , (2.22) 

where a general function of Aq is written as a Laplace-transform / (Aq) = J ds f (s) e~''^°. This 
technique has been successfully applied to sophisticated traces appearing in functional renormal- 
ization group calculations [11, 13, 21, 22], where it allowed to generalize the computations to 
operator traces which are not accessible by standard heat-kernel technique. 
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Before closing this section, it is worth noting that a non-local generalization of (2.20) appeared 
in the literature [23]. Although the result covers only to the first order in the curvatures only, it 
shows that the generalization of (2.20) to a non-local heat-kernel expansion like (1.3) is possible. 
Notably, such a generalization would give many insights to important physical problems, like the 
computation of radiative currents. 

2.4 Tensor fields on a curved background 

In this subsection, we show how to obtain the heat-kernel for the special case of a Laplace operator 
acting on scalar-, vector- and symmetric tensor fields from the general formulas (2.15), (2.16) and 
(2.20). In this case we set E = and choose = 0, such that the field strength becomes 
F^J.u = [V^, V,y] acting on scalars, vectors and symmetric tensors, respectively. 
From the commutators acting on the field spaces of interest we compute 

[V^,V,]7;" = i?^,V^: (2.23) 

where 0, v'^ and h'^^ are test scalar-, vector- and symmetric tensor fields respectively. From (2.23) 
the following relations between field strengths and Riemann curvatures are implied 

F I — 

Ml' I scalar ~ " ' 

-^Ati^°/3 1 vector ~ ^tJ.i'"' P ) (2.24) 

i? , — 2E (^'^^ 

tensor ' 

Once (2.24) are computed, we can substitute them into (2.20) and obtain the heat-kernel ex- 
pansion for the special cases. The diagonal heat-kernel expansion and its coefficients are defined 
by 

3 10 



(2.25) 



where r labels the three tensor space of interest and we adopted the basis (A.l). The coefficients 
of (2.25) are listed in table 2 together with the results of the next section. 



3 Seeley-deWitt expansion for transverse vector fields 

In this section we compute the early-time expansion of the heat-kernel resulting from a projected 
Laplace operator acting on transverse vector fields, up to third order in the curvature. To simplify 
matters, we will consider abelian vector fields setting the connection = for the remainder of 
this section. In this case = V^j and we use the former notation throughout. The final result 
(3.22) holds for compact manifolds without boundary. Throughout the computation we will make 
extensive use of the off-diagonal heat-kernel coefficients derived in the previous section. 
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3.1 Transversal vectors 

In Yang-Mills theory on flat space it is a common practice to decompose the fluctuations of 
the vector field into its transversal and longitudinal parts [24]. In Lorentz-gauge the gauge- 
transformations solely act on the transversal part, so that the latter constitutes a gauge degree 
of freedom of the theory. Following [9], the decomposition 

v^' = e + ^^(|), d^e = (3.1) 

can be implemented by introducing the projection operators 

d d'^ d d" 

which satisfy 11^ = 11 and Hh/i'^du = d^, d^Ui.fj,'^ = . In terms of these 

a^(/> = nL/t;., e = ^Tf.''v^. (3.3) 

In the context of gravity or quantum field theory on a curved spacetime, it is desirable to 
extend this decomposition to a general Riemannian (background) manifold. In this case, one has 
the local decomposition [12] 

v^' =^/^ + D^'(^ , D^i^' = . (3.4) 

This split is unique up to a constant shift in (p, which constitutes the zero- mode of D^. The 
generalization of the projectors (3.2) implementing this decomposition is straightforward 

Hl/ ^ -^M^^' ' = + D^^D-" . (3.5) 

They project v^^ to its irreducible components 

^T^'v, = , Ui^^^v, = D,,(t> . (3.6) 

In order to construct a well-defined operator trace, it is crucial to observe that the eigenvalue 
equation for the standard Laplace operator is, in general, incompatible with the transverse con- 
dition. This can be seen as follows: Suppose that is an eigenfunction of A on the space of 
transverse vector fields 

ACf = Aief, Ai>0, Z)^er = 0. (3.7) 
Applying to the left-hand-side of this equation gives 

D^Aef = [D^, A] ef = -D^^R^^^'i: (3.8) 

where we used the transverse condition together with (B.4). The divergence of the right-hand-side 
of (3.7) vanishes identically, however, so that we obtain the condition 

Di'R^,i^ = Q. (3.9) 
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On a general Riemannian manifold there is no reason that this identity holds. Thus on a Rie- 
mannian manifold the space of eigenfunctions of A does not decompose into the direct sum of 
transverse and longitudinal vector fields. 

An important special case are the Einstein manifolds discussed in appendix A. 2. These have 
the special property that R^i, oc g^^, so that (3.9) is satisfied. In other words the Einstein- 
condition ensures [A, IIt'^i/Jw'^ = which can be verified easily with the use of the equations 
(3.26). Thus, for this special case, the space of eigenfunctions of A does decompose into the 
direct sum of transverse and longitudinal vector fields. 

In order to give a well-defined meaning to the heat-kernel on the restricted set of S,^ configu- 
rations satisfying the constraint D^^^^ = 0, we define the projected Laplace operator 

A^ = nT^^AHT". . (3.10) 

The projectors ensure, that A^jy propagates the transversal modes only. Moreover, even for the 
case of a general Riemannian manifold, all eigenfunctions of A^i, are transversal vector fields 
by construction. The projectors entering into the definition (3.10) turn A'^jy into a non-local, 
pseudo-differential operator, however. Thus it is a priori unclear to which extend the standard 
results for the Seeley-deWitt expansion carry over to the case of vector fields satisfying differential 
constraints. This will be investigated in the next subsections. 

3.2 Structural aspects of the heat-trace 

Based on the previous discussion, we are now in the position to give a well-defined meaning to 
the heat-trace on the space of transversal vector fields, Trixe"*^, as the traced heat-kernel of the 
projected Laplacian A: 



'IT 



Tri Ht e-'^TAHT _ (31^^ 



The additional projector in front of the exponential serves the purpose to remove the contribution 
of the longitudinal vector modes from the zeroth order of the exponential series. Our aim is to 
compute the early-time expansion of Sit up to third order in the curvature. In this context, it 
is important to note that since A is a non-local operator and in particular not of generalized 
Laplace-type, the results of the previous section cannot be applied straightforwardly. 

The explicit computation of Sn starts from expanding the exponential appearing in (3.11) 
in a power series 

^iT = V ^—^ Tri Ht (AHt)" . (3.12) 
^-^ n! 

n>0 

Here we used the idempotency of IIt to remove the right set of projectors appearing in (IItAIIt)"'. 
Subsequently, we apply the commutation relations between the projectors and A to collect all 
powers of the (unprojected) Laplacian. Resumming the power series for each type of commutator, 
eq. (3.12) can be expanded systematically. To lowest order in the commutator this expansion has 
the form 



5iT = Tri 



s Tri 



Ht [A, Ht] e-^^ +... (3.13) 
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and the dots are symbolic for higher order commutators. Each commutator carries fixed power 
of (derivatives) of the curvature tensor. For instance, the lowest order commutator is 



A,nT/ 



(3.14) 



The higher order commutators then contain more powers of the curvature tensor or its derivatives. 
Thus eq. (3.13) constitutes a curvature expansion which allows to evaluate 5it up to a specified 
order in TZ. 

In order to apply the off-diagonal heat-kernel technique to this curvature expansion, we next 
have to deal with the inverse powers of the Laplacians appearing in the projectors and their 
commutators. Here we employ the Schwinger-trick 



1 

A" 



t 



n-l 



By combining all exponentials including the Laplace operator by applying further commutator 
relations if necessary, the traces in (3.13) can be cast into the form 



Tn 



A" 



t 



n-l 



O e-('^+*)^ 



dt. 



(3.16) 



/o ("-!)! 

where O = TZD^^ . . . D^^ constitute non-minimal operator insertions of the form (2.18). Thus the 
operator traces on the right-hand-side can be evaluated via the off-diagonal heat-kernel technique 
of section 2.3. 

The crucial new feature appearing in eq. (3.16) is the auxiliary t-integral, which can be traced 
back to the non-locality in the projection operators. Upon inserting the i/-functions (2.20), these 
integrals assume the generic form 



I{n,k-^) : 



f 



n-l 



ril-k 



(3.17) 



{n-iy.^^ ■ r(f -A;) 

The gamma functions appearing here show the typical form of an IR divergence as seen within 
dimensional regularization. Indeed the r.h.s. of (3.17) becomes singular for n + k > d/2, that 
is, if the order of the heat-kernel expansion 1Z^~^^ exceeds d/2. One would thus expect that the 
coefficients in the Seeley-deWitt expansion of Sit will develop divergences at 

0(7^'^/2)_ With 

respect to the explicit computation in the next subsection, we will use (3.17) to regularize the 
Seeley-deWitt coefficients by analytically continuing in the spacetime dimension d, assuming that 
it is sufficiently large for all integrals to converge. 

Before closing our discussion, we remark that it is necessary to keep the spacetime dimension 
unspecified, so that no terms are lost in the traced heat-equation 



Tr 



Aii-(s) 



ds 

(4vr) 



Tr 



-d/2 



[h{s)\ 



n>0 



Integrating this relation should reproduce the standard heat-kernel 



Tr[i7(i 



Tr[Ai7(s + t)]dt, 



(3.18) 



(3.19) 



which requires that none of the factors (| — Ji) vanishes in the above formula. 
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3.3 Evaluation of the heat-trace 

After discussing the structural aspects related to the evaluation of (3.11), we proceed carrying 
out the explicit computation including all terms up to third order in the curvature expansion. 
Our first task is to complete the curvature expansion (3.13) by including all (multi-)commutators 
which give rise to contributions up to order TZ'^. To find the coefficients multiplying the basis 
monomials (A.l), it is sufficient to keep track of all operator insertions which contain up to three 
powers of the curvature or terms where two covariant derivatives act on two curvature tensors. 
Terms where three or four covariant derivatives act on one curvature tensor are, in principle also 
of order 7^^, but give rise to surface-terms only and may thus be discarded. From (3.14) we 
learn that [A,nT] oc TZ. Taking an additional commutator [A, [A,nT]] then gives terms with 
two powers of the curvature or one covariant derivative acting on R/^u- From this systematics we 
conclude that the highest multi-commutator contributing to 0{TI^) is the quadruple commutator 
[A, [A, [A, [A,nT]]]]. In addition one has to keep all products of lower order commutators, that 
satisfy the above criteria. Starting from (3.12) and working out the combinatorics, the curvature 
expansion of 5it up to 0{Ti^) is given by 

SiT-E^(-^r^i^' (3-20) 

n=0 

where 
^i? =TVinTe-^^ 
5g =TVinT [A, Ht] e-*^ 
5g=TrinT ([A,nT]' + [A, [A,nT]]) e-^^ 

=TrinT (^2[A,nT][A, [A, 0x1] + [A, [ A , Ht ]] [ A , Ht ] 

+ [A, [A, [A,nT]]] + [A,nT]')e-^^ 
=TnUT (^3[A,nT][A, [A, [A,nT]]] + 3[A, [A,nT]]' 

+ [A, [A, [A,nT]]][A,nT] + [A, [A, [A, [ A , Ht]]]] ) e"^^ , 
and ~ indicates that the right-hand-side is actually a curvature expansion. 

(n) 

The actual evaluation of the S^rp can then be carried out along the lines of the previous 
subsection. Since the computation is rather lengthy and not very illuminating, we relegate these 
technical details to appendix D. The final result for Sn is obtained by substituting the interme- 
diate results (D.5), (D.8), (D.12), (D.16) and (D.18) into (3.20). In terms of the basis (A.l) it 
takes the form 



(3.21) 
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(3.22) 
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with the coefficients given in the last column of table 2. Eq. (3.22) constitutes the main result 
of this section. 

Commenting on the result for the Sit trace, first we notice that it has poles in even spacetime 
dimensions. Their origin can be traced back to the f-integrals in (3.17), that are finite for 
sufficiently large d only. We observe that the coefficients at order TZ" diverge for n > d/2. Thus, 
the poles occur in one order higher then estimated via power counting, so the highest d poles 
at each order in the curvature cancel, allowing us to analytically continue the results to lower 
dimensionality. Unfortunately, this phenomenon does not occur for all the poles and we are left 
with the bounds d > 2 and d > 4 for the orders TZ^ and TZ^. A more detailed discussion on these 
poles will be given in the next subsection. Notably, all results can be safely continued to any odd 
dimensionality. This is owed to the fact that, for odd d, the gamma functions in (3.17) never 
become singular. 

The coefficients of the transverse vector trace, given in table 2, do not factorize into the 
difference of those of the vector and the scalar trace. This is what one would naively expect 
by simply counting the degrees of freedom, but is obviously wrong because of the presence of 
"interaction terms" between transverse and longitudinal modes. Interestingly, however, there is 
a partial decoupling of the degrees of freedom at any order in the curvatures. In fact, it is always 
possible to write any coefficient of the transverse trace as the aforementioned difference modulo 
a correction: 

cP- (transverse) = (vector) — (scalar) + Qi (d) /Q2 (d) , (3.23) 

where Q2 (d) is a polynomial of one degree higher than Qi (d). This means that, in the limit 
d —7- 00, the full decoupling is realized and only the number of degrees of freedom matters, a 
situation that is reminiscent of many mean field computations at large d. 

3.4 Transversal vector fields on Einstein-spaces 

As already discussed in section 3.1, the Einstein-spaces reviewed in appendix A. 2 constitute a 
special class of Riemannian manifolds where the projection operators (3.5) commute with the 
Laplacian. In this case the basis monomials in (3.22) are not independent but related by the 
geometrical identities (A. 13). The corresponding heat-kernel coefficients are readily obtained 
from the general result (3.22) by using (A. 13) to express all curvature monomials in terms of the 
Einstein-space basis 

•51x1^^ = / d'^a:^g{id-l) + ^^^^s£' 

j_ „2 ( 5d3-7d2-58d-180 c'i i rf-16 c2\ f3 24) 

^ * \^ 36M2 "^1 180 ^^2 ^ \ ■ ) 

j_ „3 f 35d^-77d^-60M'^-3512d-75m c3 , 7d^-llld-127 c3 , 17d-269 c-3 _ d-19 c3\ \ 
* 45360d3 '^l 7560d ^^2 "i" 45360 1620 4^ J • 

Remarkably, in this case we find heat-kernel coefficients that are finite in any dimension d > 0. 
Evoking the Einstein-space limit the heat-kernel coefficients in (3.22) combine in such a way that 
all the singularities appearing for even dimensionality d cancel. 
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At this stage, it is illustrative to rederive (3.24) by using the Einstein-space condition from 
the very beginning. Since [AjHtJIes = 0, all S^j^ , n > 1 in the expansion (3.20) vanish and 

S'itIes ~ "^iT (3.25) 
is exact. Exploiting that on an Einstein-space 

Z5,/(A)0 = /(A + |)Z)„<^, Z)„/(A)^;° = /(A-|)I)«^;°, (3.26) 

the trace (D.l) can be cast into standard form without resorting to the expansion (D.2) 



oo 



+ / TVi 



e"* d Z)^Z)''e-("+*)^ 



dt . (3.27) 



Here the Ricci scalar in the exponential results from the summation of the multi-commutators in 
(D.2). The curvature expansion of the traces can then, again, be found via the off-diagonal heat- 
kernel. The additional exponential thereby leads to an exponential suppression of the integrand 
for large values of t, rendering the auxiliary integrals finite. This computation confirms the 
result (3.24). Notably, the exponential ensures that the heat-kernel coefficients appearing in 
Sixlgg are free from dimensional poles to all orders in TZ^"'\ Moreover, the fact that the general 
computation restricted to the Einstein-space limit and the Einstein-space computation lead to 
the same finite result indicates that the origin of the dimensional poles is not due to the use of 
the early-time expansion of the off-diagonal heat-kernel. We will elaborate this point further in 
the next subsection. 

We close this subsection by evaluating (3.24) on a spherical background. Using (A. 15) we 
obtain 

Sphere ~ Ws)^ J ^V9|la " Ij + gd + 360d^(d-l) ^ ^ 

I - lild^ -331d-^ -3825d-'^ -676d^ +10992d-7560 „3 nSl 

"I 45360d^(<i-l)^ ^ j ' 

which for d = 4 becomes 



(3.28) 



Sit 



Sphere,.=4 = W^j '^'^^{^ + " TSO«'^' " Jio^'^'} " (3-29) 

This trace still contains the contribution of the constant scalar mode, which does not contribute 
to the T-decomposition (3.4). In order to obtain the final result, this mode has to be removed by 
hand, leading to a modified heat-trace [25] 



R 

IT e = itiT e 

R (3.30) 



TrWe""^ =T\^iT e'^^ + e'T 



Sphere, d=4 

Here the second term is precisely (minus) the contribution of a constant scalar. Taking this 
correction term into account, the heat-kernel coefficients for transversal vectors on the sphere 
become 

c'o = 3 , 4 = i , 4 = , 4 = - 362880 • (3.31) 

Notably, this result is in complete agreement with previous computations [6, 25, 26]. 
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3.5 Regularizing the poles 

The infrared divergences of the auxihary integrals (3.17) appear in the heat-kernel expansion of 
S'lT in the form of poles in even dimensions, as shown in table 2. The first divergence in d = 4 
spacetime dimensions appears at order TZ'^, when the boundary term Ai? is neglected. Notably, 
however, the coefficients of TZg and TZ^q are finite for any dimensionality. 

Any attempt to regularize (3.22) has to be implemented with care in order not to affect the 
unambiguous result for Einstein spaces, since in this limit the poles cancel with a factor {d — 4) 
forming in the numerator. If no pole is present, these contributions would otherwise go to zero. 
In order to take the limit d — )• 4, we expand in e = d — 4 and make use of the formula 



^"-'/'^ = + /(4) - ^ log (^) ) ) + 0{e) 



n-2 



/'(4)-i/(4)logm +0(£) 



(3.32) 



where the j pole is absorbed into a (divergent) redefinition of the scale Sq = s^e^^^ . In this way 
s'q plays the role of an infrared cutoff and allows us to obtain the limit as 



+ ^'((m + Mlog(t))^? -(2sf6+^log(t))^2 

+ (3i6+T34log(t))^i -S + M5l°g(t))^4 

+ (^ + ilog(^))7^i -(|i, + ^log(f))7^3 



(3.33) 



- MJ log (f )) ^? - (lio + M log (i)) 



67 I 1 ^3 



15120 '^9 108 '^10 

The appearance of the logarithmic terms is an immediate consequence of the dimensional reg- 
ularization employed. The Seeley-deWitt expansion is still "weakly" valid after dimensional 
regularization, because the singular log(s) terms are always multiplied by a power of s, and so 
maintain a continuous s — )• limit. ^ Using this regularization scheme, in principle no information 
is lost by regularizing by a finite choice of Sq in (3.33). In fact, it is still possible to recover the 
correct Einstein-space limit from (3.33), provided that the explicit dependence of Sq on e = d — 4 
is taken into account, so that it can combine with the d-dependence of the Ricci- and scalar 
curvature tensors in their Einstein-space limit. 

To demonstrate the origin of the divergence in the low momentum contribution of the inverse 
Laplacian in the projectors (3.5), we consider a one-parameter family of "projectors" 

til, = g^,y + aD^^D^, (3.34) 



®The very same happens also in the d — s> 2 hmit. 
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that are idempotent for a = 1 only. With this definition the modified projected Laplacian can be 
expanded in the form 

Af,^ = {U^ AU^)f,^ = Ag^^ + a{2-a)Df,D^ + 0{n,A-^) , (3.35) 

where all the pseudo-differential contributions due to the commutators of the inverse Laplacian 
with the covariant derivatives are schematically contained in 0(7^, A^^). At zeroth order in 
the curvature, the heat-kernel of A is singular at a = 1, due to the degeneracy of the operator 
[2, 27]. Since for other values of a 11^ is not a projector anymore, we can exclude the degeneracy, 
that is specifically its infinite number of zero eigenvalues corresponding to the volume of the 
subspace of longitudinal modes, as source of the divergence. Tracking this modification in the 
traces (3.21) reveals a simple polynomial dependence of (3.22) on the parameter a. Therefore 
a discontinuity for the special value a = 1 is not present. An explicit calculation in fact shows 
that the dimensional poles vanish only for a = 0, reproducing the heat-kernel of the standard 
Laplacian. 

Alternatively, we can remove the zero eigenvalue from the spectrum of the Laplacian specifi- 
cally by shifting it with an infrared scale m?, defining the operator 

n;^^' = g^u + D^{A + rr?) '^D, . (3.36) 

The scale m? can be chosen arbitrarily small, suppressing the low momentum modes while leaving 
the high momentum spectrum, essentially, unaltered. Using this modified operator in place of 
the projectors in (3.11), the exponentiation via (3.15) would be done on the regulated inverse 
Laplacian (A -|- m^)~^. This accounts for an additional regularizing factor e"*'"^, which renders 
the modification of the integrals (3.17) convergent at any order of the curvature. We conclude 
that the long range (infrared) modes of the pseudo-differential projection operators cause the 
breakdown of the heat-kernel expansion (3.22). This situation is analogue to the case of a scalar 
field on curved spacetime, whose propagator can be defined via the heat-kernel in a local expansion 
only if it is massive (or otherwise IR regulated). 

For practical purposes, it is conceivable to use a modified operator like (3.36), if can be 
identified with an infrared scale already present in a particular problem. In general however, such 
a procedure has the disadvantages that the limit — t- is discontinuous and arbitrary powers 
of the dimensionless combination sm? will occur in heat-kernel traces. Instead, we suggest to 
follow (3.32) and to introduce an infrared scale in a purely dimensionally regulated setup. 

4 Summary 

In this article, we used the de Witt-algorithm to recursively determine the off-diagonal heat-kernel 
expansion of a Laplace operator on a general bundle, including an arbitrary endomorphism, up 
to third order in the curvature tensors. The algorithm has been implemented on computer 
algebra systems [17, 18], making the manual handling of the huge expressions for higher orders 
unnecessary. Our results generalize previous calculations [11, 14, 15, 16] by leaving the internal 
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Table 2: The traced heat-kernel coefficients in the early-time expansion (2.25) on the space of 
scalars, vectors, symmetric 2-tensors and transversal vectors (IT), respectively. The dimensions 
of the field spaces are denoted by di = d,d2 = ^d{d + 1) and dn = d — 1, respectively. 

space and spacetime dimension unspecified. In particular the given in eq. (2.16) is, to our 
knowledge, the most general determination of this coefficient, since it is given without tracing 
over the internal space. All the results for the traced expansion of the heat-kernel of the Laplace 
operator acting on scalar-, vector- and symmetric tensor fields are conveniently summarized in 
table 2. 

Our results allow the systematic evaluation of the traced heat-kernel, including insertions of 
covariant derivative operators inside the trace. These results become relevant for the evaluation 
of heat-traces of non-minimal differential operators, for which the heat-kernel expansion is harder 
to obtain and known to less extent [23]. In this case, an expansion and resummation of the non- 
minimal part allows to reduce the heat-trace to terms of the form (2.22) [2, 13]. This is particularly 
useful in the context of quantum field theory, where expressions of this form generally appear 
when a gauge-invariance is present. As an example, the newly computed non- minimal traces in 
section 2.2 provide the means to extent previous approximations of the functional renormalization 
group flow of gauge theories and gravity in numerous ways. In particular, they will be an essential 
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ingredient in extending the universal renormalization group machine [11] to the third order in 
the curvatures. 

As a non-trivial application of our results, we studied the traced heat-kernel on the space 
of transverse vector fields satisfying D^v^ = 0. Notably, the Laplace operator projected on 
the transverse vectors subspace becomes non-local due to the projection operators. Owed to 
this non-locality, the Seeley-deWitt coefficients of the heat-trace, on a general (i-dimensional 
manifold, diverge starting from order TZ'^^'^'^^. The divergences are of infrared nature and appear 
as dimensional poles, when dimensional regularization is adopted. Most remarkably, we showed 
that these singularities cancel, if the spacetime is an Einstein-space. The reason is that the 
Einstein condition ensures that the standard Laplacian commutes with the projection onto the 
transverse subspace. 

Remarkably, the occurrence of the singularities is not related to the use of the early-time 
expansion of the heat-kernel employed in section 3. Based on power-counting arguments, we 
estimated that divergences are present also if a non-local expansion (1.3) of the projected trace is 
employed. This estimate leads to the same divergence structure in the form of dimensional poles, 
as that indicated by the results reported in table 2. We take this as a strong indication that our 
poles are not an artifact of the expansion, but rather a genuine feature of the projected traces. 

We close our summary by commenting on the consequences of these divergences in the context 
of the functional renormalization group approach to quantum gravity. In this context, the method 
of transversal field decomposition was first advocated in [25] and subsequently employed by a 
number of groups [6, 26, 28, 29, 30, 31] to diagonalize the propagator. This technique works 
as long as the background manifold satisfies the Einstein condition, but our results indicate 
that it does not generalize to the case of generic backgrounds. Instead of using the transverse 
decomposition, one can perform a resummation of the non-minimal derivative terms appearing in 
the inverse propagator [2] . This idea has already been used to reobtain the 1-loop effective action 
of four-derivative gravity in [13], and can be generalized to a full non-perturbative computation 
of the effective action. 
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A Curvature monomials on Riemannian manifolds 

An essential ingredient for the early-time expansion of the heat-kernel is a basis for the inde- 
pendent curvature monomials at a given order 7^". In this appendix, we therefore review the 
results of [15, 32, 33] which are essential for our construction. Throughout this appendix, we will 
consider the case of d-dimensional Riemannian manifolds without boundaries. 
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A.l Curvature basis for general manifolds 

The heat-kernel expansions (2.25) and (3.22) require a basis for aU curvature invariants build 
from six or less covariant derivatives. Following [32] this basis contains 15 elements, which we 
choose as follows^ 



7^o = 


1, 










n^ = 


R, 










n = 


R\ 


nl 


— Rp,uR^ ) 




— R aRP^""^ 


nl = 


RBR, 




= R^yHiR'^^ , 




= R', 


nl = 


RR^vRP'^ , 


n 


— -Tt^j JT^ -TLq, , 




p r? J)uavB 


7^f = 






— /? E^'-'^P"! R'^ a 

— ^f-fiu^^ a/37 ' 


nl 


— R Po"E> al3n uu 

— ^^fiu -"-per ^''ap 




~ /J u^^ p a^^ a p • 











(A.l) 



Here □ = g'^'Df.D^. 

Notably, the results obtained from the off-diagonal heat-kernel techniques are, a priori, not 
given with respect to this basis, but also include additional curvature monomials, which are 
related to the ones above through the first and second Bianchi identities 

R^.[ua^] = (A.2) 

and 

A very useful collection of curvature identities implied by these identities has been given in [33] , 
and for completeness we summarize the ones important for our construction. 
At order TZ"^, (A.2) implies 

R^""'^ RliauP = ^Rp.ualBR^'^"'^ , (A. 4) 

while at TZ^ there are three relations between different contractions of the Riemann tensor 

r)p.aul3 r> o per _ 1 p pan a(3 tj uu 

Ral3^^ R"'^^'^ Rppcru = jRpu'"^ Rpcr°'^ Ra/S'^" , (A-5) 

R „WVf'^R'" /3 — _ If? per/? "'P R a^^" R" ^ RP " RP a 
^^p.avp^'' a p — 4-'i-p!/ -"-per ^^ap i p v^^ p cr-"' a p • 

In addition the combination of (A.2) and (A. 3) allows to derive 



RpuaP^R^"^^ —4:Rpu;al3R^°"^^ + '^RpuRajB-y^'^^^^ ~ Rpu'^" Rpu^^ RaP^'^ 
- AR"" 1^ RP ^ RP a 



(A.6) 



and 



RpuR^"'^ a — hR\puR^^ + Rp'^ Ru"' Ra'^ — RpyRapR^'^^^ ■ (^-7) 



^Taking total derivatives into account, there is one more invariant at 0{TZ^) and seven additional curvature 
monomials at 0(71'^), see eqs. (A. 8) and (A. 9) below. 
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Once the curvature monomials appear under the volume integral, the condition of working 
on a manifold without boundary allows us to freely integrate by parts. At 'R? this eliminates the 
surface term 

j d'^x^UR = 0, (A.8) 
while at order IZ^ we obtain the seven additional identities 

J d'^x^ R^^.^pRf^'^''^ = j d'^x^ [nl - \nl -izl + 7^3] , 



d'^x^R.^R'' = - j d'^x^lZl , (A.9) 

j d'^x^ Ra^.^R'^^'f' = - 1 d'^x^nl, 
j d'^x^R^p.^R'^f^'P = j d^x^ [_l7^3-7^3 + 7^3] , 
j d'^x^ R^,a(S;pR^''''''''' = J d'^x^ [nl - 47^i + 47^i - 47^3 - 2nl + nl + 47^?o] . 

A. 2 Curvature basis for Einstein-spaces 

A special class of Riemannian manifolds are Einstein-spaces, where the Ricci tensor is proportional 
to the metric 

R^.u = yf.uR■ (A. 10) 

In connection with the Bianchi-identities (A. 3) this definition entails 

R■,^^ = 0, iia/3;M = 0, Ra^jf,;^ = , (A.ll) 

for any dimension d ^ 2. 

As a direct consequence of these additional relations, the general basis of curvature monomials 
(A.l) degenerates, so that, up to 0{TZ^), an Einstein-space has only eight distinguished curvature 
monomials. Explicitly, these can be chosen to be 

£^ = 1, 
£^ =R, 

(A.12) 

Ci — n , C2 — J^^uaf)^ 1 

cS rj3 cS p p Tjfii/al3 c3 p po" p D /"^ pa /3 p/i 1/ pp a 

Ci — n , C2 — ixn^yaprL , C3 — n^y iXpc Ji-a/j , C4 — n ^ p a /3 ■ 

Using eqs. (A. 10) and (A.ll), the additional basis elements in (A.l) can be expressed in terms of 
theC 

'^2 — d'-'l ' 

7^f = o, 7^i = o, 'R-l = l£!, (A.is) 

'^5 ~ d^'-'l ' 'H ~ ' '^8 — d^2 • 
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Lastly, the maximally symmetric d-spheres S'^ pose a special class of Einstein-manifolds. In 
their case, the relations 

R/ii/ = ^d/J^i^R J Rfii/pa = d{d—l) i9fJ.p9iycr ~ 9lJ.(T9iyp)R ^ (-^-l^) 

determine all curvature invariants in terms of the constant Ricci scalar or, equivalently, the radius 
of the sphere. Thus, the basis of curvature monomials contains only one element at each order 
TZ"": i?". The additional basis elements contained in (A. 12) then satisfy 

c2 2 p2 

' " '^'-'^ ' (A.15) 

c3 _ 2 o3 c3 _ 4 p3 c3 _ d~2 p3 ^ ^ 

^2 - d{d-i)^ ' ^3 - W{d^^ ' - d^{d-iy^-^ ■ 

These relations are used to compare our results to previous computations of the heat-kernel on a 
spherical background in subsection 3.4. 



B Auxiliary commutator identities 

This appendix gives the explicit expressions for multi-commutators that appear in the evaluation 
of the traces over constrained vector fields in section 3. In this appendix we identify = and 
define the Laplacian A = —g^^'^D^j_Dy. Symmetrization is with unit strength (a/3) = ^ (a/3 -|- /3q). 
In order to lighten our notation we furthermore introduce the n-fold commutator 

with [ 1 ^ ] ~ -^A* ^'^'^ [ -^M ' ^ ] 1 ~ [ 1 ^ ] • Finally, we use ~ to indicate that the 
corresponding equations are exact up to terms which do not contribute to the curvature expansion 
in section 3. 

B.l Simple commutators acting on scalar and vector fields 

The computation of the operator traces (3.21) requires the explicit expressions for the multi- 
commutators (B.l) acting on scalars and vector fields up to order n = 4. Most conveniently, the 
higher-order commutators can be obtained recursively from the lower order ones. For complete- 
ness, we collect the corresponding expressions in this subsection. 

Up to order n = 4, the multi-commutators acting on scalar fields are given by 

[D^, A] ^(j) = [r^'^RJ - {AR/) + 2i?^"'^L»,) Dp ^ , 
[D^,A]^^ = (r^'^R.^RJ - {A{R^''Rj)) +2{R^^Rj),^D'^ 

fB 2) 

- {AR^'')Rj + AR^J^D^D'')Dp(t} + 2R^^'"[DpD^,A](l), 
[D^, A] ^ct) ~ 4i?^("'«^ (2R\D^DpD, - 2R^'' ^^D^D^^D, 

+ R.^^D^D^Dp - iR%\D^DfjD^')(P + 0{VJI'^) , 



23 



with the commutator appearing in the third term given by 



Notably, the first three commutators are exact, while in the fourth one we only displayed the 
terms up to 0{TZ^). 

The (contracted) commutators acting on vector fields can be calculated along the same lines. 
In order to prevent the formulas from becoming too lengthy, we thereby display terms up to 
C(7^3) only: 

+ {D,,{R^'^R"P))(t)p + (Z)„Ai?-^)<^;3 , 
[D", A] 3 (2i?„^i?"^'^'^ - R^'^R^'' - i?'^^ + (Ai2'^'^)) [R/Dpcl), - 2RP^%Dp<l^„) (B.4) 

- 2^°^^'^ [D^D^, A] 0, - m'"'-^"'D„D^D^4>, , 

A]^ ~ - AR'''''^'^(2R/DxD^D^cI), + R^^DxDpD,,(t>, - AR^\^DpD^D„(t>x 

- 2R^%^DpD^D„<j)x - m^'^^^DpD^Dxcpu - 2R^'' p^D„D^Dx<l), 
These expressions are completed by the commutator appearing in the three-fold commutator 

[D^D^, A](/), = R/D^Dp^, + R^pD^D^cP, - 2R/,'' D^Dp(l)^ - 2R^%'' DpD^cP, 

- 2R^P,''DpD^(j), - Re,p:pDPcp, + Rap;pDPcp, + Rpp;aDP<l), 

(B.5) 

- RaufDp(j)p + Ra^;uDp(l)p - Rpu-PDe^cPp + RpP,^Da(t)p 
"^Rp^ u p(j)cr R^vf a't' p ~l~ Rfi^ ■,i/a4'p • 

For further commutator relations, we refer to the appendix of [11]. 
B.2 Commutators involving composite operators 

For two continuous linear operators X, Y the Baker-Campbell-Hausdorff formula allows to write 

oo ^ 

e^ye-^ = j:-(-ir[y,X]„, (B.6) 

n=0 

where the n-fold commutator has been introduced in eq. (B.l). Identifying X = =FsA, this 
formula gives an exact curvature expansion for the commutator between the exponential of a 
Laplacian and an arbitrary operator Y 



oo -. 
n=l 



_ n! 

(B.7) 
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Using the formal properties of a Laplace-transformation, (B.7) can be extended to the com- 
mutator of Y with a function of the Laplacian 



[Y,f{A)] = [Y, A]^/W(A) 



n! 

n=l 

where f^"^^ denotes the n-fold derivative of /, and all Laplacians have been moved to the very 
right and very left, respectively. 

Based on this result, one can obtain the curvature expansion of the projection operator (3.5) 
in terms of the n-fold commutators 



nL/^. = -{f;(^(-ir[^.,A]„) D^r 



n=0 

oo 



n=0 ^ ^ 



(B.9) 



Again, this expansion has the advantage that all Laplacians have been moved to the very right 
or very left, respectively. The explicit expressions for the multi-commutators appearing in these 
expansions up to order n = 4 are given in eqs. (B.2) and (B.4). These features are essential when 
evaluating the operator traces including projectors in appendix D. 



C More derivatives of the off-diagonal heat-kernel 

In this appendix, we give some additional explicit formulas for the derivatives of the off-diagonal 
heat-kernel coefficients An{x,y) in the coincidence limit. While interesting in their own right, 
these formulas allow to cross-check computer-algorithms. 

The terms containing five and six symmetrized derivatives of Aq read 

15 1 

and 



D^sD-yDpD^D,D^)Ao 



5 5 1 g (C'-2) 

+ g-^(7/3;7-^i^M;a) + j^RiS-yR/BaRup.) + QR{SjR\p\l3\e\aR^ u p,) 
g 

+ ^-^^(<5|e|7-R l3\X\aR u\p\p) ■ 

The expressions for three symmetrized derivatives acting on Ai and one derivative acting on A2 
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contain both the endomorphism and curvature of the non-trivial bundle. They are given by 
1111 3 

Di^aDyD^j^^Ai = - -E.^^^a) + —R-i^^ya) - -^ER^^^.f^) - -E.j^aR^iu) - -^Epf^^.P^a) 

13 11 
111 1 

1 „ „ 1 



and 



D,A2 = ^i^E),, + ^E,,E + ^EE,, + ^E,pE% + ^F%E,p + ^EE^..^^ 

-I- J_ ppo" F -I- — F -I- — F F^"^ -I- — EP"^ E —F^E"^ (C 41 

^ ^ PP'^ ^ P^'P" 30 ^^''^ 45 P"'^^^ go '^''^ ^ v^-^'' 

Lp'^ /? —FPR -A- — EP"^ /? —E RP'^ A 3_ TPP m 

36^'^ 30^^M^;P+30^^ «pm;- qq^pm;-^ + 180 "'^^ 

~ -^Epa-.yR^P"'^ — —{AR)-a + -;^RR;p — Yso'^'"^'^'"^'^ ~'~ i^^^'^^'^ ^P'^^'^'P ' 

respectively. Recall that -E and F^i/ may be matrix-valued with respect to the internal bundle 
and therefore, in general, do not commute. 



The coefficients D^^pDyDaDj^-^Ai and L)(^D^)^2 which, following table 1 also enter into the 
recursive construction of ^3(2;) are not given, since they are very lengthy expressions and therefore 
of little practical value when written explicitly. If needed these coefficients may be computed from 
a supplementary Mathematica file, which will be provided on request. 



D Traces containing projection operators: intermediate results 

In this appendix, we collect the intermediate results entering the computation of the traced heat- 
kernel on the space of transversal vector fields in section 3. In this course, we will give the explicit 
expressions for the S^^ defined in (3.21) in the following subsections. 

D.l Evaluating Sf^ 

Substituting the explicit expression for IIt given in (3.5) into and employing the Schwinger- 
trick the trace can be written as 



5 



(0) 
IT 



Tn 



+ 



T\-i 



dt. 



(D.l) 



Here the first trace is a standard heat-trace on the space of unconstrained vector fields and is 
readily evaluated by substituting the heat-kernel coefficients given in the second column of table 
2 into (2.25). 



26 



The second trace is evaluated via the off-diagonal heat-kernel. We first combine the two 
exponentials using the Baker-Campbell-Hausdorff formula 



oo ^ 

L J ^-^ n\ 



n=0 



where 



(D.2) 



(D.3) 



and [D^, A]„ denotes the n-fold commutator (B.l). The commutators appearing in (D.3) ensure 
that (D.2) constitutes a curvature expansion of the trace. Following the argument at the beginning 
of subsection 3.3, all terms contributing to the basis (A.l) are generated by the first five terms 
of this expansion, so that we truncate (D.2) at n = 4. Substituting the commutators (B.2), it is 
a straightforward application of the off-diagonal heat-kernel to evaluate the T^*^) . We find 

1 



r(o) 



"(4^(s + t))'^/2 



-rpz _ a— 4 -TP 2 

2(s+t) 12 144 ''^1 ^ 360 '^"i 360 ''''3 

5^+12 ^:. 3 _ d+36 T?3 _ d+12'rj3 , d-SO-rjS , 16d+345 T?3 

140 '^l 70 '^2 108 '^3 ^ 90 '^i ^ 945 
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4d+207 ^3 d -p3 ^-174^^3 17^-102^^3 d-6 -773 
315 '^6 90 '■''7 630 ''''8 3780 '^9 + T35 ''''10 



r(2) 

^(3) 
2.(4) 



+ 



(47r(s + t))'^/2 
s + t 

1a' 



2{s+t) • 



3/V2 

6l?3 I 8x'3 I lx'3 I 19x'3 22't->3 i 56x'3 i 1 -pS 4 t-jS 
5 '"<-l "T 5 ''^2 "T 6 ''*'3 "T T5 '"*^4 ~ 15 '"^b "T T5 ''*^6 "T 15 ''*^7 ~ T5 ''^5 



(D.4) 



(47r(s + 1))'^/2 



2( 



i_nl + ^ (37^? -nl + 2nl - (m\ 



(4vr)'i/2(s + t)rf/2+i 



(4vr)"!/2(s + t)d/2+2 

Based on these results, the auxiliary t-integration in (D.l) can be executed. Adding the vector- 
trace contribution, we arrive at the final result of the form 



S 



(0) 

IT 



1 



(47rs)rf/2 



10 



(D.5) 



i=l i=l 

with the (i-dependent coefficients given in the second column of table 3. 
D.2 Evaluating Sf^ 

We proceed with the evaluation of s'^r^ defined in (3.21). Substituting the commutator (3.14) 
and using the orthogonality of the projectors IIt • Hl = 0, this trace becomes 



S- 



(1) 

IT 



-TVi 



(D.6) 
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S- 



w 

IT 



A 

72 



d-1 

i d+6 

1 d+10~ 
72 (i-2 



d I d^-32d+180 
180 180d(d-2) 



d-1 
180 



12 



1 5d^+32d+464 
^ 120 1680{d-4){d+2) 

J d^+40d^-64d-1120 

30 



5- 



W 

IT 





1 

d(d+2) 
1 

Jd+T) 





d^+4d^+24ii-24 
'6(d-2)d((i+2)((i+4) 

d^+4d-24 
6(d-2)d(d+4) 

1 

6(d-2)(i 
d^+d+lO 



5 



IT 



1 



d(d+2) 
1 





-10 



336 



d 
840 



1296 



+ 



840(d-4)d(d+2) 

d+14 
1296(d-4) 

d2-30d+236 



2835 ^ 30 ^ 
945 " 



1080 ' 1080(d-4)((i-2) 

16d-*+377d^ - 1954d^ +2272d-30240 



J_ 

30 



11340(<i-4)(d-2)d(d+2) 

4d^+223d-1460 
3780{d-4){d+2) 



d 


1 


d+2 


1080 


72 


1080{d-4) 


d 


1 


d-1 72 


7560 


90 


7560{d-4) 




17(d-l) 
45360 


1 
180 




d-1 


^ 90 




1620 



6(d-2)d(d+2) 

__df+8d+32 
6d(d+2)(d+4) 

d^-12d-32 
3(d-2)d(d+2)(d+4) 









d"'-14d+8 
2(d-2)d(d+2)(d+4) 

4 

(d-2)(d+2)(d+4) 

d^ + 12d-4 
6{d-2)d{d+2)(d+4) 

d'^+14d^-4d+40 
6(d-2)d(d+2)(d+4) 

d^-2d+16 
3d(d+2)(d+4) 

4(d^+8) 
3{d-2)d{d+2){d+4) 









(n) 

Table 3: Traced heat-kernel coefficients of the partial traces S\r^ entering the computation in 
section 3.3. 



In order to cast this expression into standard form, we express the projection operators via (B.9), 
where the inverse Laplacians already appear to the very left and very right of the commutator 
insertions. Since (D.6) already contains one explicit power of the curvature, the series of higher- 
order commutator can be terminated at order n = 3: 



c(i) ^ 
•^IT 



I 



n=0 / 



-sA 



. n=0 



(D.7) 

Collecting the inverse powers of A and employing the Schwinger-trick, s'^^ assumes the form 
(3.16). Substituting the explicit expressions for the commutators (B.2) and (B.4), the trace 
becomes accessible by the off-diagonal heat-kernel. Evaluating all subtraces with Mathematica, 
the final result for 5^^^ is of the form 



?(i) 
'it 



(47rs)'^/2 



10 



4 = 1 



(D. 
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with the expansion coefficients given in the third column of table 3. One can verify that 

vanishes on an Einstein-space, since [A, IIt] = for this special case. We find that the coefficients 

remain finite in d = 4, so that they do not contribute to the leading divergences present in 
c(0) 

D.3 Evaluating S^^ 

(2) 

The 5|rp-trace (3.21) naturally decomposes into two subtraces containing the product of two 
single commutators (3.14) and the double-commutator [A, [AjIIt]]. The latter can conveniently 
be constructed recursively from (3.14). The operator insertions in these traces are then of the 
schematic form A~"C'A~"'CA^"". In order to collect all Laplace operators in a single function 
one has to take into account that the curvature tensors are not covariantly constant and thus 
do not commute with the (inverse) Laplacians. These terms can be taken into account via the 
commutation relation 



(D.9) 



which captures all terms contributing to the required order. 

With this formula, the explicit results are found via Mathematica and read 



TViHt [A, li-rf e 



(47rs)rf/2 
+ 



d{d+2) I'M 



dizl) 



^3-4^^+32^+40-1? 3 ^ d3+8d^-4d-32 y^3 _ 3^3+2^^-24^-32 -7:. 3 



{d-2)d{d+2) 



6(<i+4) 



2(d+4) 



3{d+4) 



3(<i+4) 



(D.IO) 



and 



TriHT [A, [A, Ht]] e"*^ 



(47rs)<'/2 



d(d+2) 



dnl) 



%[d-2)d{d+2) 



3(<i3+20d-16) ^3 12(d-4)(d+2) ^3 , q/^ , o^^:.3 
2+4 '"^-l 2+4 ''^2 + ^l" + ^J'"*^3 



2{d^ + d+ 10)7^| + 



d+i 

3 _L 4(d-2)(d2+4d+16) .^3 _ 2(d-4)(3d^ + 10d+16) ^3 



d+4 



d+4 



(D.ll) 



respectively. 

(2^ 

The final result for S\r^ is the sum of (D.IO) and (D.ll) and has the curvature expansion 



^(2) 
'it 



1 



(47rs)'^/2 



d'^x^ 



10 



4 = 1 



1=1 



(D.12) 



with the coefficients given in the fourth column of table 3. 
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D.4 Evaluating and Sl^ 

The contributions of and can be computed along the hues of the previous subsections. 
In order to present the final result, it is convenient to introduce the abbreviations 

1 1 



.(4) 



Ci 



(47rs)'^/2 d{d + 2){d + 4:) 



d'^x^ d{d - i)n\ + [d^ - s)nl - 2nl 

+ Mill + d{d + 4)7^| - 2d{d + 2)7^^ 



(D.13) 



and 



C2 = 



1 



d'^x^ {d - 6)nl + 2{d + 2)nl + 2dnl - 2dTZ, 



(D.14) 



(47rs)'^/2 (d + 2)(d + 4) 

Notably both of these combinations vanish identically on an Einstein-space. The subtraces ap- 
pearing in S^^ are given by 



(D.15) 



2Tri 


Ut 


[A, 


Ut] 


[A, 


[A 


Ut]] 


e" 


sA 


= 2Ci, 


Tri 


Ht 


[A, 


[A, 


Ut] 


[A 


, Ht] 


e" 


sA 


= -Ci 




TVi 


Ut 


[A, 


A, 


A, 


Ht]]] 


e" 


sA 


= -C2 








Tri 


Ut 


[A, 




e~ 


sA 


= 0. 



Substituting these results into (3.21) yields 



S- 



(3) 
IT 



C1-C2. 



(D.16) 



For the traces contributing to we find 



3 Tri 



Tn 



Ut [A, nT][A, [A, [A, Ht]]] e 



3Tri 



Ut [A, [A, Ht]]' e 



Ut [A, [A, [A, nT]]][A, Ht] e" 



Tn 



Ut [A, [A, [A, [A,nT]]]] e 



-sA 



-sA 



■ -3s~^C2 
■2s-^C2. 



(D.17) 



The final result for is the sum of these subtraces and reads 

S^r^ = 3 S ^ C2 • 



(D.18) 
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